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ABSTRACT 

Using observations from the Chandra and XMM-Newton X-ray observatories, we examine the 
interaction between the intra-group medium and central radio source in the nearby NGC 4261 
galaxy group. We confirm the presence of cavities associated with the radio lobes and estimate 
their enthalpy to be ~2.4x 10 58 erg. The mechanical power output of the jets is ^10 43 erg s -1 , 
at least a factor of 60 greater than the cooling luminosity in the region the lobes inhabit. We 
identify rims of compressed gas enclosing the lobes, but find no statistically significant tem- 
perature difference between them and their surroundings, suggesting that the lobe expansion 
velocity is approximately sonic (M < 1.05). The apparent pressure of the radio lobes, based 
on the synchrotron minimum energy density argument, is a factor of 5 lower than that of the 
intra-group medium. Pressure balance could be achieved if entrainment of thermal gas pro- 
vided additional non-radiating particles in the lobe plasma, but the energy required to heat 
these particles would be ~20 per cent, of the mechanical energy output of the radio source. 
NGC 4261 has a relatively compact cool core, which should probably be categorised as a 
galactic corona. The corona is capable of fuelling the active nucleus for considerably longer 
than the inferred source lifetime, but can be only inefficiently heated by the AGN or con- 
duction. The expansion of the radio lobes has affected the structure of the gas in the galaxy, 
compressing and moving the material of the corona without causing significant shock heating, 
and expelling gas from the immediate neighbourhood of the jets. We discuss the possible im- 
plications of this environment for the duration of the AGN outburst, and consider mechanisms 
which might lead to the cessation of nuclear activity. 

Key words: galaxies: individual (3C 270, NGC 4261) — intergalactic medium — galaxies: 
active — cooling flows — X-rays: galaxies 



1 INTRODUCTION 

The importance of radio galaxies as sources of heating in groups 
and clusters of galaxies has become increasingly apparent in re- 
cent years. Active galactic nuclei (AGN) are now considered the 
most likely mechanism acting to balance radiative cooling of the 
hot intra-group medium (IGM, McNamara & Nulsen 2007; Peter- 
son & Fabian 2006). Approximately 50 per cent, of galaxy clusters 
(e.g., Sanderson et al. 2006) and perhaps as many as 85 per cent, of 
galaxy groups (Dong et al. 2010) have a cooling region centred on 
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a giant elliptical or cD galaxy, and the link between cool cores and 
nuclear activity in these galaxies is well established (e.g., Burns 
1990; Mittal et al. 2009; Sun 2009). 

X-ray and radio imaging has provided a means to examine the 
interaction between radio galaxies and their environment, and re- 
vealed that many nearby groups and clusters contain complex struc- 
tures associated with the radio jets and lobes. Shocks and sound 
waves driven by the expansion of radio jets may directly heat the 
gas (e.g., Nulsen et al. 2005; Fabian et al. 2006; Sanders & Fabian 
2007; Forman et al. 2007) while cavities may inject energy by do- 
ing work on the gas, or reduce its ability to cool by lifting material 
out from the densest, most highly enriched regions (e.g., Fabian 
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et al. 2005; Wise et al. 2007; Blanton et al. 2009). The enthalpy of 
cavities associated with radio lobes has been shown to be sufficient 
to balance radiative cooling in many systems (Birzan et al. 2008), 
and there is evidence that in some systems shocks may deposit sim- 
ilar amounts of energy into the gas (e.g., Million et al. 2010). 

However, the mechanisms by which AGN are fuelled, and 
thus the processes involved in starting and stopping nuclear activity 
are not clear. Correlations between the AGN power output and the 
Bondi accretion rate suggest that Fanaroff-Riley class I (FR I) radio 
sources could be fuelled by hot gas (Allen et al. 2006; Hardcastle 
et al. 2007; Balmaverde et al. 2008). While actual Bondi accretion 
would require unrealistic efficiencies (McNamara et al. 201 1), this 
at least suggests that fuelling by gas cooled from the hot phase is 
a possibility. Fuelling by large reservoirs of cold molecular gas or 
neutral hydrogen appears to be less likely, with poor correlation 
between cold gas mass and radio power (McNamara et al. 2011; 
Emonts et al. 2007). To form a feedback system, the AGN must be 
able to effectively heat its fuel supply, located in its immediate sur- 
roundings. Shocks originating in the nuclear region could be effec- 
tive in this regard, while cavity heating would require conduction 
to distribute the energy azimuthally and perhaps radially. 

While these mechanisms may be feasible in many systems, 
they appear difficult to achieve in those dominant ellipticals which 
host galactic coronae, cool cores with radii of only a few kilopar- 
secs. Coronae are thought to be largely unaffected by the radio jets 
of the AGN they host, and conduction into the cool gas from the 
surrounding IGM is suppressed (Sun et al. 2007). Owing to their 
small size, relatively small heating efficiencies (<1 per cent, in 
some cases) would raise the temperature of the corona gas to that 
of the ambient medium, suggesting that conduction is strongly sup- 
pressed by magnetic fields, and that the jets tunnel through the cool 
core with little or no interaction. Coronae are thought to consist of 
gas lost from stars within the cool region. Their cooling rate is of- 
ten similar to the rate of stellar mass loss in the core (Sun et al. 
2007), implying that if the AGN is fuelled by cooling gas, a corona 
can fuel its AGN over long periods without requiring inflowing 
gas from the IGM. Destruction of coronae by mergers or by ram- 
pressure stripping seems infeasible given the relatively large num- 
bers of such systems observed (Sun et al. 2007). Coronae thus raise 
serious problems for the feedback model, since they can apparently 
provide enough gas to fuel an AGN for long periods, without be- 
ing heated by the jets this activity produces. While this possibility 
is as yet speculative, there is at least one example of a corona fu- 
elling an unusually old radio galaxy (O'Sullivan et al. 2010). AGN 
in corona systems are also known to have different radio properties; 
the radio power of radio galaxies in large cool cores scales with the 
luminosity of the cool core, whereas those in coronae show no such 
relationship (Sun 2009). 

The question of how the supply of gas to the AGN can be 
stopped is important in all cluster and group-central AGN, and 
particularly so in coronae. Addressing this issue is difficult in that 
it requires examination of the AGN/gas interaction on very small 
scales. Many cavity systems studied to date are distant, making 
such an examination difficult, and this is particularly true of corona 
systems. It would therefore be useful to study nearby radio galaxies 
embedded in small-scale cool cores or coronae, to determine what 
effect the radio source has on the core, and how this might affect 
fuelling of the AGN. 

One such system is the NGC 4261 group, which has a well- 
established cool core and whose central elliptical hosts the well 
known FRI radio source 3C 270. Our aim in this paper is to use 
the available deep Chandra and XMM-Newton observations of the 



group to examine the interaction between the radio source, the cool 
core and the surrounding intra-group medium. We can thus deter- 
mine the effects of the ongoing AGN outburst (or series of out- 
bursts) on the gas, the energies and timescales involved, and the 
current physical structure of the gas. NGC 4261 is particularly 
suited to such a study, since the well-defined axis of its jets re- 
moves any question of large projection effects in determining the 
size of structures, and its radio properties and black hole mass are 
relatively well constrained. We describe the group, galaxy and radio 
source in more details in Section 1.1. In Section 2 we describe the 
Chandra and XMM-Newton observations and their reduction, and 
in Section 3 the results of imaging and spectral analysis of these 
data, including estimates of the age and expansion rate of the cav- 
ities. We further analyse these results in Section 4 to examine the 
energy output of the AGN, the properties and particle content of the 
radio lobes, and the state of the cool core. These results and their 
relationship with, and implications for, other group- and cluster- 
central radio galaxies are discussed in Section 5, and we summarise 
our conclusions in Section 6. 

Throughout this paper we assume Ho=70 km s _1 Mpc -1 and 
a distance for NGC 4261 of 31.6 Mpc, in line with that adopted by 
Worrall et al. (2010). This gives an angular scale of 1' = 9.2 kpc. We 
note that Humphrey et al. (2009) use a smaller distance (29.3 Mpc) 
based on surface brightness fluctuation measurements (Tonry et al. 
2001). Our conclusions would not be altered were we to adopt 
this distance value. Worrall et al. (2010) also adopt the redshift 
of NGC 4261, 0.00746 (Trager et al. 2000), as appropriate for 
their goal of studying the AGN and jets. As our aims concern the 
IGM, we adopt the mean redshift of the group, 0.00706 (Nolthenius 
1993). However, we note that adopting the galaxy redshift would 
not significantly alter the fit results. 

1.1 NGC 4261 

NGC 4261 and its surrounding group have been extensively stud- 
ied, and it is useful to review what is known of the system before 
moving on to discuss our analysis. NGC 4261 is a boxy, cuspy- 
cored , slowly rotating E2 galaxy (Ravindranath et al. 2001), whose 
primary axis of rotation is close to the major axis, leading to the 
conclusion that the galaxy is prolate (Davies & Birkinshaw 1986). 
The galaxy hosts a low-ionisation nuclear emission-line region 
(LINER, Ho et al. 1997) and the bright FRI radio source 3C 270 
(~19 Jy at 1.4 GHz, Kiihr et al. 1981), whose twin jets lie close 
to the plane of the sky. The galaxy has a kiloparsec-scale kine- 
matically decoupled core (Cappellari et al. 2007) which contains a 
100 pc-scale disk of dust and cool molecular and atomic gas (Jaffe 
et al. 1993; Jaffe & McNamara 1994) whose rotation axis is closely 
matched with the axis of the radio jets (Ferrarese et al. 1996). On 
large scales, deep optical imaging has revealed a faint tidal tail to 
the northwest and tidal fan extending southeast from the galaxy 
(Tal et al. 2009) and there is evidence of anisotropy in the globular 
cluster distribution (Giordano et al. 2005). These disturbed features 
suggest that NGC 4261 underwent tidal interactions or merger with 
another galaxy within the past 1-2 Gyr. 

Diffuse X-ray emission was first detected in the NGC 4261 
group using the ROSAT PSPC, which revealed a gaseous halo ex- 
tending to at least 40' (>360 kpc) with a luminosity weighted mean 
temperature of ~0.85 keV (Davis et al. 1995; Worrall & Birkin- 
shaw 1994). ASCA observations confirmed the presence of a cen- 
tral, spectrally hard X-ray source (Matsumoto et al. 2001). The 
poor spatial resolution of these data prevented any detailed exami- 
nation of structure within the halo, but more recent observations by 
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Chandra and XMM-Newton have revealed X-ray features which 
are clearly related to the radio source and its interaction with its 
environment. On scales of a few arcminutes, XMM-Newton imag- 
ing showed arm-like features enclosing the inner edges of the radio 
lobes (Croston et al. 2005), and confirmed that at least in the west- 
ern lobe, these correspond to the edges of a cavity (Croston et al. 
2008). 

On smaller scales, a short Chandra exposure revealed X-ray 
jets corresponding to the inner few kiloparsecs of the radio jets 
(Gliozzi et al. 2003; Zezas et al. 2005), strongly indicating that the 
jet axis is close to the plane of the sky. Analysis of a more recent, 
longer observation suggests a synchrotron origin for the jet X-ray 
emission, and finds wedge-like regions of reduced surface bright- 
ness along the jet axis (Worrall et al. 2010). These are interpreted 
as conical volumes in which the thermal gas has been displaced by 
relativistic plasma from the radio lobes, suggesting that up to 20 per 
cent, of the gas within ~10 kpc of the galaxy core has been moved 
by the action of the AGN jets and lobes. Radial spectral analysis 
confirms that the group temperature profile follows a typical form, 
with a central cool core (kT ~0.6 keV) of radius ~10 kpc, at the 
upper end of the range observed in other corona systems, and a tem- 
perature peak at ~200" (kT ~1.6 keV), with a decline at larger 
radii (Helsdon & Ponman 2000; Humphrey et al. 2006). These pro- 
files have been used to place limits on the total mass of the group 
(~6x 10 13 M Q , Humphrey et al. 2006), and on the central super- 
massive black hole (~4.4x 10 8 M ) which are in good agreement 
with dynamical estimates (Humphrey et al. 2009; Ferrarese et al. 
1996). 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Chandra 

NGC 4261 has been observed twice by the Chandra ACIS instru- 
ment, first during Cycle 1 on 2001 May 26 (Obsld 894), for ~35 ks, 
using 1/2 subarray mode, and again in Cycle 9 on 2008 February 12 
(Obsld 9569) for just over 100 ks. A summary of the Chandra mis- 
sion and instrumentation can be found in Weisskopf et al. (2002). In 
both observations the S3 CCD was placed at the focus of the tele- 
scope and the instrument operated in faint mode for the first obser- 
vation and very faint mode for the second. We have reduced the data 
from the pointings using CIAO 4.1.2 and CALDB 4.1.3 following 
techniques similar to those described in O' Sullivan et al. (2007) and 
the Chandra analysis threads 1 . The level 1 events files were repro- 
cessed, very faint mode filtering was applied to the second dataset, 
bad pixels and events with ASCA grades 1 , 5 and 7 were removed, 
and the cosmic ray afterglow correction was applied. The data were 
corrected to the appropriate gain map, the standard time-dependent 
gain and charge-transfer inefficiency (CTI) corrections were made, 
and a background light curve was produced. The observations did 
not suffer from significant background flaring, and the final cleaned 
exposure times were 29.6 and 100.9 ks. While data from the entire 
detector were examined, for the purposes of this study we primar- 
ily use the S3 CCD, as the radio source and the galaxy fall on that 
chip. 

Identification of point sources on S3 was performed using the 
CIAO task WAVDETECT, with a detection threshold of 10~ 6 , chosen 
to ensure that the task detects < 1 false source in the field, working 
from a 0.3-7.0 keV image and exposure map from the combined 

1 http://asc.harvard.edu/ciao/threads/index.html 



observations. Source ellipses were generated with axes of length 4 
times the standard deviation of each source distribution. These were 
then used to exclude sources from most spectral fits. An extended 
source was detected coincident with the peak of the diffuse X-ray 
emission; this was not excluded. 

Spectra were extracted using the SPECEXTRACT task. Spec- 
tral fitting was performed in XSPEC 12.6.0. Abundances were 
measured relative to the abundance ratios of Grevesse & Sauval 
(1998). A galactic hydrogen column of 1.75xl0 20 cm~ 2 , drawn 
from the [ftools] task nh and based on the survey of Kalberla et al. 
(2005), was adopted in all fits. This differs slightly from the col- 
umn of 1.58X10 20 cm~ 2 adopted by Worrall et al. (2010), based 
on the survey of Dickey & Lockman (1990). The difference prob- 
ably arises from the finer angular resolution of the Kalberla et al. 
survey (0.675° compared to 1°), but testing suggests it has no sig- 
nificant effect on our spectral fitting. 90 per cent errors are reported 
for all fitted values. Spectra were grouped to 20 counts per bin, and 
counts at energies above 7 keV and below 0.7 keV (see below) were 
ignored during fitting. 

Background spectra were drawn from the standard set of CTI- 
corrected ACIS blank sky background events files in the Chandra 
CALDB. The exposure time of each background events file was 
altered to produce the same 9.5-12.0 keV count rate as that in 
the target observation. A region enclosing the AGN and jets was 
excluded from the estimation of the count rate of the target ob- 
servation, to avoid any contamination from source photons. Very 
faint mode background screening was applied to the background 
data sets where appropriate. Comparison of source and background 
spectra suggested some mismatch between the source and back- 
ground spectra, mainly below 0.5 keV. This is not unexpected, as 
the soft X-ray background arises largely from hot gas in our galaxy, 
and from coronal emission associated with solar wind interactions, 
and thus is both spatially and temporally variable (e.g., Kuntz & 
Snowden 2000; Snowden et al. 2004). There are also indications 
that the spectral shape of the background has changed since the cre- 
ation of the blank-sky background files (c.f. the ACIS background 
cookbook ), which could contribute to the disagreement at low en- 
ergies. NGC 4261 also lies on the outskirts of the Virgo cluster 
and close to the galactic north polar spur (Bohringer et al. 1994). 
It is possible that emission from both these sources could contami- 
nate our observations. We therefore ignored energies below 0.7 keV 
when performing spectral fitting, so as to avoid any biases arising 
from inaccuracies in estimating the soft X-ray background. 

2.2 XMM-Newton 

The NGC 4261 group was observed with XMM-Newton during Cy- 
cle 1 (2001 December 16) for just over 33 ks (Obsld 0056340101) 
and again in Cycle 6 (2007 December 16 and 18) for a total of 
130 ks (Obslds 0502120101 and 0502120201). We reduced all 
three observations, but as Obsld 0502120101 has by far the longest 
exposure (~127 ks before filtering) our analysis focussed on this 
dataset. 

The EPIC instruments were operated in full frame mode, 
with the medium optical blocking filter. A detailed summary of 
the XMM-Newton mission and instrumentation can be found in 
Jansen et al. (2001, and references therein). Reduction and anal- 
ysis were performed using techniques similar to those described in 
O'Sullivan et al. (2007). The raw data from the EPIC instruments 

2 http://asc.harvard.edu/contrib/maxim/acisbg/COOKBOOK 
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were processed with the XMM-Newton Science Analysis System 
(SAS V.9.0.0), using the EPCHAIN and EMCHAIN tasks. Bad pix- 
els and columns were identified and removed, and the events lists 
filtered to include only those events with FLAG = and patterns 
0-12 (for the MOS cameras) or 0-4 (for the PN). The total count 
rate for the field revealed significant background flaring. Times 
when the total count rate deviated from the mean by more than 
3cr were therefore excluded. The effective exposure times for the 
MOS and PN cameras were 72.9 and 46.8 ksec respectively for 
Obsld 0502120101. 

Images and spectra were extracted from the cleaned events 
lists with the SAS task EVSELECT. Response files were generated 
using the SAS tasks RMFGEN and ARFGEN. The central AGN is rel- 
atively X-ray bright, and a significant number of out-of-time (OOT) 
events are found on the PN detector, visible in the PN data as a 
trail extending from the centre of the source toward the CCD read- 
out. An OOT events list was created using EPCHAIN, and scaled 
by 0.063 to allow statistical subtraction of the OOT events from 
spectra and images (c.f. the XMM-Newton users handbook' ). Point 
sources were identified using EDETECTXHAIN, and regions corre- 
sponding to the 85 per cent encircled energy radius of each source 
were excluded. A source corresponding to the active nucleus was 
not excluded. 

Creation of background images and spectra for the system was 
hampered by the fact that the group X-ray halo extends beyond the 
field of view. Use of the "double-subtraction" technique (Arnaud 
et al. 2002; Pratt et al. 2001) involves correcting blank-sky expo- 
sures to match a source-free area of the observation; for NGC 4261 
this is not feasible. We therefore adopt several approaches. For de- 
projected radial spectral profiles, we use a local background ex- 
tracted at the extreme edge of the field of view and discard results 
from the outer annulus of the radial profile as potentially biased. 
The effect of the deprojection should act to reduce any bias in the 
successive inner annuli. Comparison with our deprojected Chan- 
dra profile and with that of Humphrey et al. (2009) shows good 
agreement, suggesting that our results are not seriously affected. 
For spectral analysis of particular X-ray structures we use local 
background spectra extracted close to the region of interest, so as 
to subtract off any overlying group and background emission. 



2.3 Very Large Array 

We analysed 1.5 GHz data in C-array configuration retrieved from 
the Very Large Array (VLA) public archive (project AL693). 
The observations were made in two 25-MHz bands centred on 
1365 MHz and 1646 MHz, in May 2008 for a total integration time 
on source of approximately 38 minutes. We used the NRAO Astro- 
nomical Image Processing System (AIPS) package for the data re- 
duction and analysis. Data calibration and imaging were carried out 
following the standard procedure (Fourier Transform, Clean and 
Restore). Phase-only self-calibration was applied to remove resid- 
ual phase variations and improve the quality of the image. The final 
image has an angular resolution of 16" x 15" and a rms noise level 
(lcr) of 0.3 mjy beam -1 . 



' 5 http://xmm.esac.esa.int/external/xmm_user_support/documentation/ 



3 RESULTS 

3.1 IGM structures associated with the radio source 

We initially examined images extracted from both Chandra and 
XMM-Newton, smoothed at a range of scales, to determine the dis- 
tribution of diffuse emission. The AGN is the brightest source in the 
field, and the jets are clearly visible in the Chandra data. The dif- 
fuse emission is strongly centrally peaked, and fills both the Chan- 
dra and XMM-Newton fields of view. 

On scales of a few arcminutes, similar to the size of the stel- 
lar component of the galaxy (see Figure la), it is clear that mod- 
erately luminous diffuse emission extends north and south of the 
central AGN, to a distance of at least 2' (~18 kpc). This extension 
is perpendicular to the jet axis, and poorly aligned with the major 
axis of the galaxy. On larger scales, fairly uniform emission ex- 
tends throughout the field of view of both instruments. However, 
comparing archival VLA 1.5 GHz radio maps to the X-ray images, 
there are hints of structures associated with the radio lobes. In par- 
ticular, after subtraction of point sources from the Chandra image 
and refilling of the resulting holes using the CIAO dmf ilth task, 
heavy smoothing reveals an apparent arc of X-ray emission extend- 
ing along the southern edge of the eastern radio lobe, and hints of 
a similar structure south of the western lobe (see Figure la). These 
could be parts of a shell of compressed gas surrounding cavities ex- 
cavated by the expanding radio lobes, but the limited field of view 
of the S3 CCD prevents us from determining whether the structures 
completely enclose the lobes. 

The XMM-Newton images are large enough to cover the 
whole radio source, but are more severely affected by the bright 
AGN emission, owing to the broader XMM-Newton point spread 
function. Other point sources in the galaxy, and sources associated 
with other galaxies in the field of view, are also problematic, and 
unresolved sources increase the level of noise in the area of interest. 
We therefore adopted two approaches aimed at clearly determining 
whether there are cavities corresponding to the radio lobes, and the 
overall form of any structures in the IGM. As a simple test, we 
applied the wvdecomp wavelet decomposition and smoothing al- 
gorithm (Vikhlinin et al. 1998) to examine the amount of structure 
on different spatial scales. Removal of features with small smooth- 
ing scales (sJ3 pixels or 13.2") allowed us to effectively subtract 
point sources, small noise features and much of the AGN emission. 
Adopting a 5a detection threshold, we obtained the image shown 
in Figure lb. This clearly shows the north-south bar of emission 
across the galaxy core, but also shows a looped structure to the east, 
and corresponding arms of emission extending west from the north- 
ern and southern ends of the bar. These correspond to the western 
cavity and small eastern surface brightness decrement identified by 
Croston et al. (2008). 

An alternate method for revealing any X-ray structures asso- 
ciated with the radio lobes is to model the AGN and large scale 
group halo, subtract these models from the image and examine any 
residual features. Such models need not be physically meaning- 
ful, so long as they provide a good approximation to the emission 
components we wish to remove. We fitted 2-dimensional surface 
brightness models using the CIAO SHERPA package (Freeman et al. 
2001). Point sources were removed from the image (as described in 
section 2.2), and a larger region was used to exclude the bright X- 
ray point source RX J1219. 8+0545. Models were convolved with 
the monoenergetic exposure map and PSF, determined for an en- 
ergy of 1 keV, approximating the mean photon energy of the data. 
A background image consisting of the scaled particle-only image 
and scaled pn OOT events image was used, and convolved and un- 
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Figure 1. a) Chandra 0.3-3 keV image, binned by a factor of 2 and smoothed with a Gaussian of sigma 16 pixels (15.8")- VLA 1.5 GHz contours are overlaid, 
starting 3<r above the r.m.s. noise level of 0.3 mJy beam -1 and increasing by factors of 2. The restoring beam size is 16" X 15.1" (HPBW). Annular regions 
used to examine surface brightness across the cavity rim are marked in white, b) XMM-Newton 0.3-3 keV wavelet smoothed MOS+pn image, with scales 
chosen to exclude point sources and emphasise diffuse emission. Dark linear features show chip gaps and bad columns, and should be ignored, c) XMM- 
Newton residual map, showing the image after removal and filling of point sources, and subtraction of the best fitting surface brightness model for the galaxy 
and larger group halo, smoothed with a Gaussian of sigma 4 pixels (17.6"). Solid ellipses indicate the optical D25 contours of NGC 4261 and several smaller 
galaxies in the field. Dashed ellipses indicate the regions used to define cavity size and position, d) XMM-Newton residual map, as in the lower left panel, 
overlaid with VLA contours and the partial annuli used to extract surface brightness profiles across the cavity rims. 



convolved flat models used to account for the X-ray background 
and any residual particle background respectively. 

Experimentation showed that modelling the AGN as a point 
source (a delta function or narrow gaussian convolved with the 
PSF) was not effective, probably because of emission from the jets 
and dense gas in the centre of the galaxy. We therefore modelled 
the core with a /3-model. The model was fixed to be circular, since 
the extended emission to north and south tended to drive the model 
to extreme ellipticities. The best fitting parameters for this model 
were r core =1.0i];'o" an d /9=0.55i2'.o2i ^ ut we a g am emphasize 
that this model is not intended to be physically meaningful, but 
simply to provide an approximation to the core emission, to allow 
it to be subtracted. A second /3-model was added to model large- 
scale group emission, with slope and core radius fixed at the values 
determined from ROSAT (r core =1.68\ ,0=0.31 Davis et al. 1995), 



but a low normalisation was found for this component, suggesting 
that the central /3-model accounted for most of the emission in the 
field. Figure 2 shows the XMM-Newton radial surface brightness 
profile, with fitted model components for comparison. 

To examine any residual structure, we subtracted the best fit- 
ting model from an image in which the point sources had been 
removed and replaced (using the dmfilth task) and the OOT 
readout streak subtracted. The resulting residual map is shown in 
Figures lc and d. The emission north and south of the core pro- 
duces the strongest residuals, but structures associated with the ra- 
dio lobes are clearly visible. In the west, arcs of emission extend 
along the south side of the radio lobe, and along part of the north 
side. There is an apparent gap, or weakening of this emission in 
the northwest quadrant. On the east, emission corresponding to the 
lobe boundary is more diffuse, but extends approximately to the far 
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Figure 2. Upper panel: 0.3-3 keV Chandra surface brightness profile of 
NGC 4261, compared to estimates of the contribution from low-mass X-ray 
binaries (LMXBs) and stellar sources (coronally active binaries and cata- 
clysmic variables) in the galaxy. The observed surface brightness is marked 
by a solid line with error bars indicating lcr uncertainties. The expected 
surface brightness from unresolved LMXBs and stellar sources is based on 
the V-band surface brightness profile of the galaxy, scaled using the rela- 
tions of Kim & Fabbiano (2004) and Sazonov et al. (2006). The shaded 
region indicates lcr uncertainties on the sum of these two emission sources. 
The dotted line indicates a Chandra 1 keV point spread function, arbitrarily 
normalised, for comparison. The 0.3-3 keV emission is clearly steeper than 
the profile expected from stellar sources and X-ray binaries (and therefore 
the optical light profile), but broader than the profile expected for a point 
source. Lower panel: XMM-Newton exposure corrected 0.3-3 keV surface 
brightness profile, compared to the fitted model used to search for resid- 
ual structure. Although a five component model was fitted, it is clear that 
the data are primarily described by a central point source, a fairly compact 
/3-model and a flat background. 



end of the radio lobe, and encloses it to north and south. Agree- 
ment with the wavelet smoothed image is generally good, though 
the eastern cavity rim is less clearly defined and appears to extend 
further east. Immediately outside these structures is a band of neg- 
ative residuals, visible as darker regions in Figure lc/d. This is an 
artefact of the surface brightness fitting, caused by the model over- 
estimating the emission at this radius as a consequence of the north- 
south bar and cavity rims. 

To test the significance of these structures we extracted 0.3- 
3 keV surface brightness profiles across the eastern and western 
ends of the lobes. We used partial circular annuli of width 5 pix- 
els (22") with limiting angles chosen to avoid regions which seem 
to be affected by unsubtracted point sources. The annuli were cen- 
tred at R.A. 12 h 19' n 14 s , Dec. +05°49'53" (J2000) for the west- 
ern and R.A. 12 h 19 m 30 s , Dec. +05°49'24" for the eastern lobe. 
The regions, overlaid on the surface brightness residual image to 
make their relationship to the rims clear, are shown in Figure Id 
and the resulting profiles in Figure 3. The difference between the 
east and west regions is clear; the western profiles have much lower 
central surface brightness than the eastern profiles. However, both 
the southeastern and southwestern profiles show a strong peak at 



Figure 3. XMM-Newton surface brightness profiles across the the outer 
quadrants of the radio lobes, using bins of width 22" (5 EPIC-pn pixels) 
as shown in Fig. Id. The eastern quadrants have higher central values, con- 
sistent with the residual image, and both southern quadrants show a clear 
peak at ~ 100", consistent with cavity rims. The northwest quadrant, where 
a break is observed in the residual image, shows no evidence of such a rim. 
The northeast quadrant has a gradual decline from a central peak, probably 
because both point sources and chip gaps affect the area where the rim is 
expected. 

~100", followed by a decline to ~200", after which the outer 
points rise again. The peak at ~100"corresponds to the structures 
at the edge of the radio lobes, which are likely cavity rims. The 
outer rise appears to be caused by point sources or small clumps 
of emission with no clear large-scale structure, or by a large resid- 
ual to the southeast surrounding RX J1219. 8+0545. We therefore 
consider the difference in surface brightness between the highest 
and lowest values, and find that the cavity rims are significant at la 
(SW) and 5.7cr (SE) confidence. 

The northwestern quadrant shows no significant change in sur- 
face brightness across the edge of the radio lobe, confirming the 
presence of the break in the rim seen in this area in the residual 
image. The northeast profile shows a slow decline from high val- 
ues inside the lobe, with no clear rim. This is probably due to the 
diffuseness of the emission in this region; the residual image also 
shows no clear rim, but a broad area of emission coincident with 
the radio and declining towards the end of the lobe. We therefore 
conclude that the lobes of 3C 270 have formed cavities in the IGM, 
and that they are surrounded by rims of gas swept up by the lobe ex- 
pansion. However, the western rim is either weak or incomplete in 
its northwest quadrant, and the complexity of the emission around 
the eastern lobe suggests that either the cavity has a complex mor- 
phology or that there are additional emission components present, 
such as unresolved point sources. 



3.2 Structure in the galaxy core 

On smaller scales, we can examine the structure of the diffuse emis- 
sion immediately surrounding the AGN and jets. Figure 4 shows a 
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Figure 4. 0.3-2 keV Chandra image of the core of NGC 4261, with 0.493" 
pixels smoothed with a Gaussian of sigma 2 pixels (upper panel) and adap- 
tively smoothed Chandra temperature map {lower panel). The colour bar 
indicates the approximate gas temperature in keV. The two panels have the 
same alignment and angular scale. 



Chandra image of the core and jets. As noted by Worrall et al. 
(2010), there are wedge-like surface brightness decrements in the 
area immediately north and south of the western X-ray jet, and a 
suggestion of similar structure on the eastern side (see also Figs. 2 
and 4 of Worrall et al. 2010). The emission bordering these decre- 
ments is somewhat brighter and extends to larger radius than the 
emission along the north-south axis, forming an X-shaped struc- 
ture. The extension to the southwest is particularly notable, while 
that to the northwest appears weakest. Worrall et al. (2010) suggest 
that the decrements are caused by expulsion of the thermal gas by 
the expanding radio lobes. In this case the X-structure would likely 
be this expelled gas, compressed and driven north and south from 
the jet axis. 

As a simple test of the importance of the different sources of 
emission in the core, we compared the 0.3-3 keV Chandra surface 
brightness profile with that expected from discrete sources in the 
stellar population of NGC 4261 (see Figure 2). We estimated the 
contribution from the low mass X-ray binary (LMXB) population 
using a y-band optical surface brightness profile (determined from 
observations described in Bonfini et al. 201 1) scaled to the expected 
X-ray luminosity of the LMXB population, based on the relation of 



Kim & Fabbiano (2004) and taking a total A" -band luminosity of 
2.2x 10 11 L K © (Ellis & O'Sullivan 2006, scaled to our chosen dis- 
tance). We note that this is likely to be an overestimate, since we 
have already excluded all resolved point sources but we expect the 
profile to be largely unchanged. The luminosity is converted to an 
expected number of counts in our band using a powerlaw model 
with r=1.8. A similar scaling is used to determine the expected 
contribution from coronally active binaries and cataclysmic vari- 
ables, based on the relation of Sazonov et al. (2006) and scaling 
using a 0.5 keV, solar abundance APEC model. Surface brightness 
profiles were determined using elliptical regions chosen to match 
the optical light distribution, with axis ratio 1.24 and position an- 
gle 68.3° (where the position angle is defined as the angle between 
the major axis and due west). Regions corresponding to the X-ray 
jets were excluded, but testing shows that these make only a minor 
contribution in the central bins of the profile. It is clear that the ob- 
served surface brightness profile is steeper than that expected for 
emission from the stellar population, but more extended than the 
Chandra 1 keV point spread function, confirming that much of the 
emission in this band arises from hot gas. 

3.3 Spectral properties of the cavity rims 

Given the ongoing jet activity in the system, it is possible that the 
radio lobes are still expanding. If they are expanding subsonically, 
their X-ray bright rims must consist of compressed gas which may 
be expected to be cooling more rapidly than its surroundings. If 
they are expanding at or above the sound speed, the lobes will 
cause shock heating of the surrounding gas. Using XMM-Newton, 
we found no significant temperature differences between the cavity 
rims, and regions immediately inside and outside the cavities. This 
is unsurprising, since the regions used are necessarily large and in- 
clude unresolved point sources and gas at a range of temperatures. 

While the field of view of the ACIS S3 CCD limits our abil- 
ity to examine the edges of the radio lobes with Chandra, we can 
examine the southern parts. Figure la shows a series of partial el- 
liptical annuli which were used to search for a surface brightness 
jump across the southeast lobe edge. These regions were chosen 
to match the shape of the apparent surface brightness feature, with 
large widths to increase the signal-to-noise ratio in each bin. Fig- 
ure 5 shows the exposure corrected surface brightness in each bin. 
A peak is found at the position of the lobe rim, where the surface 
brightness exceeds that immediately outside the rim at >3a signif- 
icance (a 3.4a difference between bins 3 and 4). This corresponds 
reasonably well with the feature observed in the XMM-Newton im- 
ages. 

If this surface brightness increase is the result of a shock, the 
standard Rankine-Hugoniot jump conditions (e.g., Landau & Lif- 
shitz 1959) can be used to estimate the shock velocity. The change 
in surface brightness of a factor ~ 1 . 1 5 indicates a shock with Mach 
number M — 1.05 ± 0.02, suggesting that if a shock is present, 
the lobe is expanding only marginally supersonically. The presence 
of a shock can only be confirmed by the detection of a tempera- 
ture increase. Spectra were extracted across the rim, using larger 
regions inside and outside to maximise the signal-to-noise ratio. 
The spectra contained ~ 1900-5400 counts, with the background 
contributing 24-35 per cent. Figure 5 shows the temperatures mea- 
sured from an absorbed APEC model fitted to these spectra, and the 
fit parameters are shown in Table 1 . The uncertainties are large, as 
expected for relatively small regions with extensive foreground and 
background group emission, and while the best fitting temperature 
of the rim is higher than that of its surroundings, the temperature 
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Radius from lobe centre (arcsec) 

Figure 5. 0.3-3 keV Surface brightness and temperature from elliptical an- 
nuli extending across the southern edge of the eastern radio lobe. Radii are 
measured along the minor axis of the ellipses, which is aligned roughly 
south. The rim of the lobe falls in bin 3 of the surface brightness profile. 
Surface brightness uncertainties are la, temperature uncertainties 90% 



difference is not statistically significant. However, the temperature 
increase expected for a M = 1.05 shock is only a factor 1.05, and 
the data are consistent with heating at this level. 

Chandra observations of the southwest radio lobe of Centau- 
rus A have shown that the X-ray emission at the boundary of the 
lobe is best characterised as a single component absorbed pow- 
erlaw model, and arises in large part from synchrotron emission 
(Croston et al. 2009). To test the possibility of a synchrotron com- 
ponent contributing to the emission around the lobes of NGC 4261, 
'powerlaw' and 'apec+powerlaw' models were fitted to the spectra 
extracted across the lobe rim. Fitted model parameters are given 
in Table 1. An absorbed powerlaw model was a poor fit to all 
spectra. The addition of a powerlaw to the absorbed apec mod- 
els did not significantly improve the quality of fits, and in all 
cases the photon index of the powerlaw was poorly constrained 
and the normalisation consistent with zero. Fixing the photon in- 
dex at r=2 (compatible with the data and similar to the values 
found for Centaurus A), we are able to place a 90 per cent, upper 
limit on the synchrotron flux from this portion of the lobe rim of 
Fx, 0.3-7.0 ^3.18x 10 -14 erg s _1 cm~ 2 . Scaling to the total area 
of the lobe rims, this is equivalent to an upper limit on luminos- 
ity of Lx,o. 3-7.0 <3.3xl0 40 erg s _1 . The equivalent (detected) 
thermal luminosity of the rims is Lx,o. 3-7.0=4.15 x 10 40 erg s _1 . 
A similar analysis of the north-south bar also found no evidence of 
non-thermal emission. 

The Chandra profiles extend across the cavity rim in a direc- 
tion roughly perpendicular to the jet axis. We might expect a higher 
expansion rate, and therefore greater surface brightness differences 
(and potentially temperature differences), at the jet tips. The lack 
of a detected rim corresponding to the northwest quadrant of the 
western radio lobe (see Section 3.1) is informative in this regard, 
and we do not see strong surface brightness features corresponding 
to the east or west extremes of either lobe. However, if the expan- 



sion were very subsonic, gas swept up by the lobes would be able to 
flow out into regions of lower pressure, and a rim feature would be 
weak or nonexistent. It is therefore likely that the lobe expansion 
rate is comparable to the external sound speed. We conclude that 
the cavity rims consist of material swept up and compressed by the 
expansion of the lobes, and at most only weakly heated by shocks 
or compression. 

3.4 Spectral map 

To examine the spatial variation of temperature in the gas, we pre- 
pared a temperature map using the technique developed by David 
et al. (2009) which takes advantage of the close correlation between 
the strength of lines in the Fe-L complex and gas temperature in 
~1 keV plasma. Most of the emission from such gas arises from 
the L-shell lines from Fe-XIX (Ne-like) to Fe-XXIV (He-like). 
For CCD resolution spectra, these lines are blended to form a single 
broad peak between approximately 0.7 and 1 .2 keV. The centroid or 
mean photon energy of this peak increases with the temperature of 
the gas as the dominant ionisation state of Fe shifts from Fe XIX in 
0.5 keV gas to Fe XXIV in 1 .2 keV gas, with Fe XVII and Fe XVIII 
providing the strongest line emission at the temperatures seen in the 
core of NGC 4261. The mean photon energy of the blended L-shell 
lines is independent of energy above ~1.2 keV. 

We can thus estimate the temperature distribution of the gas by 
mapping the mean photon energy in the 0.7- 1 .2 keV band. To deter- 
mine the relation between kT and mean photon energy we require 
information on the detector response and the typical properties of 
the gas. We therefore extracted a spectrum from a ~1' region cen- 
tred on the galaxy, excluding the central point source and jets, and 
fitted it using an absorbed apec model. Based on this fit, we chose a 
value of 0.5 Zqsls representative of the abundance and temperature 
of the central part of the group, and simulated spectra based on this 
model, with redshift set to that of NGC 4261 and the hydrogen col- 
umn set to the galactic value. These showed that for temperatures 
between ~0.5 and ~1.1 keV, the relationship is approximately lin- 
ear (kT — —4.43 + 5.68 < E >). Prior testing in other systems 
suggests that the relationship is relatively insensitive to variations 
in abundance and hydrogen column (O'Sullivan et al. 2009). 

Previous studies have shown the core temperature of 
NGC 4261 to be 0.6-0.7 keV (e.g., Worrall et al. 2010) and the IGM 
temperature to be ~ 1.6 keV (e.g., Humphrey et al. 2009). The core 
is therefore suitable for mapping using this technique. The lower 
panel of Figure 4 shows the resulting map. The estimated tempera- 
tures are comparable with those derived from previous spectral fits 
(and our own radial temperature profile, see Section 3.5). 

We do not expect to see the jet in the temperature map, since 
its powerlaw emission is spectrally flat compared to the strongly 
peaked plasma component, and indeed no features corresponding to 
the jet are observed in the map. The X-structure observed in surface 
brightness is clearly visible, which confirms that the core gas dis- 
tribution is disturbed and that the arms of the X consist of cool gas 
with a temperature ~0.65 keV. The regions of low surface bright- 
ness along the west jet have a higher temperature (0.8-0.9 keV), in 
agreement with the temperature measured by Worrall et al. (2010). 
On slightly larger scales gas with temperatures <0.95 keV appears 
most extended to the northeast and southwest. The cavity rims and 
the north-south bar are not observed in the temperature map, con- 
firming that these do not contain significant quantities of cool gas. 
We conclude that NGC 4261 hosts a small cool core, whose struc- 
ture has been disturbed by the expansion of the radio lobes and 
which is at present not spherically symmetric. 
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Table 1. Fit parameters for the Chandra spectra extracted across the southern sector of the east lobe. 



Model 



Parameter 



Bin 



1 



wabs*apec 



wabs*(apec+powerlaw) 
with T=2 



wabs*powerlaw 



kT (keV) 
Abundance ( Zq) 
X 2 I d.o.f. 

kT (keV) 
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X 2 I d.o.f. 

r 

X 2 / d.o.f. 
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1.36 



+0.21 
-0.07 



50+ ' 35 
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1 35+°' 19 
^ -0.07 
C3+O.35 
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2.15 1 



;+0.19 
-0.18 
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The powerlaw photon index was fixed at T=2 in the apec+powerlaw model, as it was poorly constrained when freely fitted. 



3.5 Radial Spectral Analysis 

In order to examine the underlying structure of the intra-group 
medium and estimate the energy required to form the cavities, we 
extract spectral profiles from the Chandra and XMM-Newton data. 
Regions corresponding to the cavities are excluded, and the core 
and jet are excluded from the Chandra spectra using a rectangu- 
lar region of length 30" and width 4" centred on the nucleus. For 
the Chandra data, annular regions were selected to contain 3000 net 
counts, and for XMM-Newton a total of 10000 net counts across the 
three EPIC cameras. The only exception is the innermost bin of the 
XMM-Newton profile, whose radius was fixed at 35" so as to sim- 
plify comparison with the Chandra profile. The spectra were fitted 
with a deprojected, absorbed APEC model. An additional power- 
law component was included to account for emission from the AGN 
in the innermost bin, and from unresolved point sources in the cen- 
tral three XMM-Newton bins. Abundances were tied between bins 
where necessary to stabilise the deprojection. Figure 6 shows the 
resulting temperature and density profiles, and the pressure profile 
derived from them. Pressure was calculated as P = nkT where 
n = 2n e . 

The Chandra and XMM-Newton profiles agree reasonably 
well at large and small radii, but there is some disagreement in 
temperature in the ~35-65" bin, and in density between ~35" and 
170". The temperature disagreement is likely caused by the differ- 
ence in the radii of the outermost bins of the two profiles, and the 
relatively low resolution of the profiles in comparison to the tem- 
perature gradient, particularly in the XMM-Newton profile, where 
the full range of temperatures across the gradient is represented by 
only a single model temperature in the deprojection. Scattering of 
low energy emission outward from the central bin is also a possi- 
bility, given the broad point spread function of XMM and the high 
surface brightness of the core. However, the large spectral extrac- 
tion region of the central bin should reduce this effect. The density 
difference may be a product of the difference in extraction regions, 
since the XMM-Newton annuli cover a larger area north of the core 
than the ACIS spectra, which reach the edge of the S3 CCD. The 
outermost annuli of both profiles are excluded since they contain 
emission from the outer parts of the group halo (which extends 
beyond the field of view) and would thus give incorrect density 
estimates. However, the level of agreement is sufficient for our pur- 
poses. We note that the Chandra profiles also agree well with those 
of Humphrey et al. (2009). 

The shapes of the temperature and density profiles are typi- 
cal of cool-core galaxy groups. Defining the cooling region as that 
volume within which the temperature declines (~10 kpc radius), 
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Figure 6. Deprojected temperature, density and pressure profiles, derived 
from the ACIS S3 (errorbars) and XMM-Newton (grey rectangles) observa- 
tions. 



the cool core is relatively small. The core has a relatively constant 
temperature of ~0.6 keV, in agreement with the estimate in the 
temperature map. The temperature peaks at ~ 1.6 keV, between 10 
and 40 kpc from the AGN, and flattens or declines at larger radii. 
The pressure profile appears relatively smooth, though this may be 
misleading since the annular average conceals azimuthal variations 
(c.f., Figure 4). 



4 ANALYSIS 

4.1 Outburst energy budget and timescale 

To estimate the mechanical energy output of the AGN, we approx- 
imate the area of the cavities with ellipses, and assume rotational 
symmetry about their major axis. The ellipses are shown in Fig- 
ure lc. Neither the shape of the radio lobes nor the shape of the 
cavities in the X-ray is straightforward, and therefore there is sig- 
nificant uncertainty in the following estimates. The east cavity in 
particular is poorly defined, and could be smaller than we have as- 
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sumed. The lobes may also be overpressured and still expanding, 
which would lead to an underestimate of the energy they contain. 

Assuming the lobes to be completely filled by relativistic 
plasma, we can estimate their enthalpy to be APV where V 
is the lobe volume and P the external pressure. The two el- 
lipses are defined as having semi-major and -minor axes of 
18.3x16.7 kpc (West) and 20.7x15.3 kpc (East). Their size is 
such that the external thermal pressure of the IGM varies from ~3- 
80x 10~ 12 erg cm~ 3 along their length. We use pressure at the mid- 
point of the lobes, from the Chandra profile, ~5 x 10~ 12 erg cm -3 , 
but note that the XMM-Newton pressure estimate is a factor of 2 
higher. We therefore estimate the enthalpy of the two cavities to be 
~2.4xl0 58 erg. 

The timescale over which cavities have formed, and therefore 
of the associated AGN outburst, is often estimated from dynam- 
ical arguments (e.g., Churazov et al. 2001; Dunn et al. 2005). In 
this system, the lobes appear still to be connected to the core, and 
their shape suggests that the jets are still the dominant factor deter- 
mining the lobe size and location, rather than buoyant forces. We 
can estimate the timescale of formation based on the time taken 
for the lobes to expand to their current size at the external sound 
speed, but as it is likely that the lobes expanded more rapidly in 
the past, and are only now approaching pressure balance with the 
surrounding IGM, this should be considered as an upper limit on 
the timescale of formation. For fcT=1.5 keV, the sound speed is 
585 km s _1 . The sound speed in the cool central region is consider- 
ably lower (380 km s _1 ) but it seems likely that expansion through 
this region will have happened at the onset of the outburst, when 
the source is likely to have been strongly overpressured. The tips 
of the lobes are ~38 kpc (West) and 45 kpc (East) from the cen- 
tral point source, suggesting that the lobes must be at most 75 Myr 
old. The mechanical power output of the AGN (including both jets) 
is therefore ^10 43 erg s _1 , with greater power output indicated if 
the radio lobes expanded supersonically over some portion of their 
history. This value is consistent with the energy output estimated 
for NGC 4261 by Croston et al. (2008). It exceeds the estimate of 
Cavagnolo et al. (2010) by an order of magnitude, but Cavagnolo 
et al. did not identify the main cavities, considering only the struc- 
tures associated with the jets in the core region. Lower limits to 
the total radiative power output of the AGN can be taken from 
the intrinsic X-ray luminosity of the nucleus (1.4xl0 41 ergs -1 
Worrall et al. 2010), or the extended radio luminosity at 178 MHz 
(!/U=10 40 erg s^or L„=24.73 W Hz" 1 Balmaverde et al. 2006; 
Chiaberge et al. 1999). 

The longer eastern lobe extends to ~5' (~45 kpc) from the 
nucleus. The total X-ray luminosity of the gas within this radius 
is ~1.6xl0 41 erg s _1 (0. 3-7.0 keV), of which approximately two 
thirds arises in the central 10 kpc. The mechanical power output of 
the AGN exceeds this value by at least a factor of 60. If the forma- 
tion timescale of the cavities is longer than we have estimated, this 
discrepancy will be reduced, but given the evidence for near-sonic 
expansion of the lobes, it is unlikely that the radiative losses can 
balance the energy output of the radio source. Even if we consider 
only the work done in moving the thermal gas as the lobes expand, 
PV, the energy involved would certainly be sufficient to balance 
cooling. 

4.2 Energy balance in the cool core 

We consider the state of the cool core of NGC 4261 by examining 
the energy balance between radiative cooling and heating. Heat- 
ing could be supplied by the AGN, by supernovae in the core, or 



through conduction from the hotter surrounding IGM. We define 
the core as including the inner six bins of the Chandra deprojected 
temperature profile, giving an outer radius of ~9.8 kpc (~1'), but 
note that Worrall et al. (2010) have estimated that up to ~20 per 
cent, of the gas in this volume may have been evacuated by the ra- 
dio source. As discussed in Section 4.1, the mechanical power of 
the jets is >10 43 erg s _1 , while the total X-ray luminosity of the 
gas in this region is ~10 41 erg s _1 (0.3-7.0 keV). The jets therefore 
need only heat the gas with 1 per cent, efficiency to balance cool- 
ing. Taking a maximum temperature of 1.6 keV for the IGM and 
a typical core temperature of 0.6 keV, we can estimate the energy 
that would be required to heat the cool core to the temperature of 
the IGM, ~9.1 x 10 56 erg. Given the timescale over which the AGN 
has been active (75 Myr at most), heating with an efficiency of 3.5 
per cent would have destroyed the cool core. In fact, the heating 
efficiency is likely to be considerably lower, since the core still has 
fcT=0.6 keV and the mechanical power of the jets will be greater if 
the timescale of the outburst is shorter than our estimated limit. 

Given the temperature gradient between the core and sur- 
rounding IGM, conduction could potentially be a significant heat 
source. Following the methods described in O' Sullivan et al. 
(2007), and working from the properties of the IGM gas imme- 
diately outside the core (bin 7 of the Chandra temperature and 
density profiles), we estimate the mean free path of electrons (in 
the absence of significant magnetic field) to be ~8 kpc. Since this 
is comparable to the scale of the temperature gradient, conduction 
into the cool core is likely to be saturated (Cowie & McKee 1977). 
The saturated heat flux per unit area into the core will be: 

/ 2kT \ 1 / 2 

q 3 at=rn e ( I n e kT erg s _1 cirT 2 . (1) 

\7rm e / 

Assuming the core to be spherical, the saturated conduction 
rate would thus be 1.8xl0 42 ergs -1 . Given the highly aspheri- 
cal structure of the core, the surface area over which conduction 
can occur will be higher, leading to an increased rate. However, 
this estimate does not account for the effect of magnetic fields and 
should therefore be considered an upper limit. While a field align- 
ment which increases conduction is possible, it is unlikely and the 
common assumption of a tangled field would reduce conduction 
significantly. We discuss this further below. 

While a factor of ~6 less than the mechanical energy avail- 
able from the jets, the conduction rate still exceeds the radiative 
cooling rate by a factor of ~18, suggesting that conduction must 
be strongly suppressed for the gas to have retained its cool temper- 
ature. Heating at this conduction rate would destroy the cool core in 
16 Myr, a shorter timescale than that required to form the cavities. 
Assuming a timescale equal to that of the AGN outburst, we find 
that conduction must be suppressed below the rate estimated from 
equation 1 by a factor of >5. 

However, this timescale is almost certainly too short. There are 
several reasons to believe that the temperature structure of the IGM 
prior to the AGN outburst must have been similar to that currently 
observed. Firstly, most nearby galaxy groups have a temperature 
structure similar to that of NGC 4261, with a central cool core, peak 
in temperature at moderate radii and temperature decline at large 
radii (Rasmussen & Ponman 2007; Sun et al. 2009), indicating that 
this is probably a moderately stable state. Secondly, the formation 
of the observed temperature gradient over a timescale comparable 
to that of the outburst would be difficult. If the AGN formed the 
temperature gradient by heating the ICM at 10-70 kpc, the energy 
required would be ~7.5 x 10 58 erg (assuming an initial temperature 
equal to that of the core, ~0.6 keV). This exceeds the enthalpy of 
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the cavities by a factor of 3, and would require strong additional 
heating of the gas, presumably by shocks, for which there is no evi- 
dence. Since the virial temperature of the group would be ~0.6 keV 
in this scenario, such strong heating would be likely to cause con- 
vective motions in the IGM, and might be sufficient to unbind gas 
from the group gravitational potential. The heated region also ex- 
tends well outside the tips of the radio lobes (70 kpc compared 
to ~45 kpc), but the timescale over which heated gas could dis- 
perse outward is longer than the AGN outburst timescale. The sonic 
timescale for gas in the hot region corresponding to a distance of 
60 kpc is ~100 Myr, and convective motions would take several 
times longer to transport the gas outward. It is therefore likely that 
the appropriate timescale over which to consider conduction into 
the cool core is much longer than that of the AGN outburst, and the 
suppression factor for conduction is correspondingly higher. 

A suppressed conduction rate is not unexpected in the cool 
core of NGC 4261. The simplest assumption of a tangled mag- 
netic field would lead to significant suppression. Magnetic field 
lines aligned radially along a temperature gradient, which would 
enhance conduction, are likely to be realigned by convective gas 
motions, taking on a tangential alignment which will strongly sup- 
press conduction. This heat-flux-driven buoyancy instability (HBI, 
Parrish & Quataert 2008) can reduce the conduction rate by large 
factors. The timescale for saturation of the instability is long and 
dependent on the steepness of the temperature gradient (Parrish 
et al. 2009), and any field alignment can be disrupted if turbulent 
motions are introduced into the gas (Parrish et al. 2010). However, 
HBI provides a potential mechanism for suppressing conduction to 
the level required. 

The number of supernovae occurring within the core can be 
estimated based on the stellar population in the region. NGC 4261 
was taken to be a prolate ellipsoid whose major axis lies in the 
plane of the sky, and which is described by a Sersic (1968) model 
with effective radius 56.81", axis ratio 0.80, and shape parameter 
n=5.44 (determined from observations described in Bonfini et al. 
201 1). These physical parameters are consistent with the limits de- 
termined from stellar velocity measurements, ignoring the prob- 
able small inclination of the galaxy away from the plane of the 
sky (Davies & Birkinshaw 1986). The Sersic model was depro- 
jected using the three dimensional luminosity density distribution 
described by Chakrabarty & Jackson (2009), normalising to the to- 
tal optical luminosity Lb=5x 1O 1o Ls0. 

Assuming a supernova rate of 0.166 per 100 yr per 1O 1O L_S0 
(Cappellaro et al. 1999), and an energy release of 10 51 erg, we esti- 
mate the heating rate within the core to be 6.4x 10 40 erg s _1 , if the 
energy release is efficiently coupled to the gas. This agrees within 
a factor of 2 with the X-ray luminosity of the core. However, since 
the X-ray surface brightness profile is considerably more peaked 
than the optical light profile (see Fig. 2), it seems likely that super- 
nova heating will be most effective in the outskirts of the cooling 
region, where it may balance radiative cooling, but that it will be 
unable to prevent cooling in the high density gas near the AGN. 

To determine whether the cool core could potentially fuel 
the AGN, we estimate the isobaric cooling time and mass depo- 
sition rate in the innermost Chandra spectral bin (radius ~0.3- 
0.8 kpc). We find a cooling time of ~60 Myr, and a mass depo- 
sition rate M C o O ;=0.033 M yr _1 . Comparing this with our es- 
timate of the mechanical power output of the AGN, we find that 
the efficiency of the AGN in converting cooling gas into energy 
would be e conv = P mech /Mc 2 =\.Sx 10~ 3 . Zezas et al. (2005) 
argued for an accretion rate much lower than this mass deposition 
rate (4.3 x 10 _6 r; _1 M Q yr -1 , based on the Eddington ratio of the 



AGN) but it is unclear whether the core X-ray emission component 
used to make this estimate arises from radiatively inefficient accre- 
tion, or from the radio jets (Zezas et al. 2005; Worrall et al. 2010). 
Based on the rate of stellar mass loss from AGB stars (Athey et al. 
2002), we estimate that ~0.095 Mq yr -1 of gas is injected into the 
core from stars. Supernovae will make only a minor mass contri- 
bution, based on the rate of 2 x 10~ 3 SN yr _1 , but will contribute 
to heating the gas to the observed temperature of ~0.6 keV. The 
total mass injection rate is similar (within the large uncertainties) 
to the rate of loss through cooling, suggesting that there is no need 
for a flow of cooling gas into the core to sustain its current size and 
gas content. It therefore seems possible for the cool core to provide 
cooling gas to fuel the AGN over long timescales, replenished only 
by gas lost from the stellar population within the core. 

In summary, it seems likely that neither the AGN jets nor con- 
duction from the surrounding IGM is effective in heating the cool 
core. Supernova heating may balance radiative cooling in the outer 
part of the core, but is unlikely to do so in the denser inner regions. 
However, stellar mass losses within the core should be sufficient to 
replenish any gas lost, allowing the core to maintain its current size 
and to continue to supply cooling gas to fuel AGN activity over 
timescales longer than that of the current AGN outburst. 

4.3 Pressure balance in the lobes 

As discussed in Section 3.3, while the lobes may be overpres- 
sured with respect to their environment, the difference in pres- 
sures is probably not large. Given approximate pressure equilib- 
rium, a comparison of the apparent pressures of the thermal plasma 
of the IGM and the relativistic plasma of the radio lobes can pro- 
vide insight into the particle content of the lobes (e.g., Hardcastle 
& Worrall 2000; Dunn et al. 2005; Birzan et al. 2008). The lobes 
of FR I radio sources are generally found to have apparent syn- 
chrotron pressures well below the IGM pressure (e.g., Feretti et al. 
1992), and this is generally assumed to indicate either that the lobes 
contain a significant population of non-radiating particles, or that 
the assumption of equipartition is incorrect, with the magnetic field 
providing additional pressure support. 

We estimate the physical properties of the radio lobes from 
the 1.5 GHz VLA archival data. The volume of the lobes is derived 
from the ellipse regions in Figure 1, as described in Section 4.1. 
We adopt an electron spectral index of p=2.2, equivalent to a radio 
spectral index of a=0.6, which compares well with the flux density 
measurements compiled in NED 4 , and the value of q=0.62±0.17 
estimated by Kiihr et al. (1981). We make the usual assumptions of 
minimum energy conditions, in which the total energy content of 
the relativistic particles and magnetic field are roughly equal, and 
a low energy cutoff in the electron energy distribution at a Lorentz 
factor jmin = 10. We can then estimate the strength of the mini- 
mum energy magnetic field B min and the pressure of the relativis- 
tic plasma P m in.en., which is defined to be 

p _ -Bmtn , (1 + k)jg e 

i mtn.en. — ~ ~r oT/^ ' ^ ' 

2/ioe oV4> 

where V is the lobe volume, E e is the energy of the electron popu- 
lation and no is the permeability of free space {An x 10~ 7 in S.I. or 
47T in cgs units). <j) is the filling factor of the lobes (assumed to be 1) 

4 The NASA/IPAC Extragalactic Database (NED) is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. 
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Table 2. Magnetic field, pressure and particle content of the radio lobes 

Property Lobe 

West East 
1.5 GHz Flux density (Jy) 8^7 SU 

Volume (cm 3 ) 6.3xl0 68 6.0xl0 68 



mm. en 



( erg cm d ) 



6.2 

io- 12 



6.2 

io- 12 



and k the ratio of energy in non-radiating particles to the energy in 
relativistic electrons (assumed to be 0). The estimated parameters, 
assuming (l+k)/0=l, are shown in Table 2. The factor e is depen- 
dent on the ordering of the magnetic field, with e=3 for a tangled 
field, and e=l for a uniform or stretched field (Leahy 1991). We 
have adopted e=3, but we note that e=l is assumed in some com- 
monly used formulae (e.g. Burns et al. 1979; Fabian et al. 2002). 
For comparison, adoption of e=l would approximately double our 
estimate of lobe pressure. 

Comparing with the IGM pressure at the lobe midpoint (~ 
5 x 10~ 12 erg cm~ 3 ) the difference in IGM and lobe pressures of 
a factor 5 agrees closely with the estimate of Croston et al. (2008). 
Under the assumption that P m i„. e n. is an underestimate and the 
lobes and IGM are in fact in pressure equilibrium, we can rear- 
range Equation 2 to determine the changes in filling factor or non- 
radiating particle component necessary to explain the pressure im- 
balance. However, we must account for the dependence of B min 
on k and (f>, and the dependence of E e on B min . In simplified 
form, B min oc [(1 + k)/>] 1/(a+3) , and E e oc £T^ +1) (Wor- 
rall & Birkinshaw 2006). We find that pressure balance is achieved 
if (l+k)/0=5.1. If, as seems likely, the lobe is still mildly overpres- 
sured with respect to the IGM, this value can be considered to be a 
lower limit. 

A reduced filling factor could occur in a number of ways; the 
lobes may contain a radio-quiet plasma component such as gas en- 
trained by the radio jets, there may be variation in the relative den- 
sity of magnetic field and relativistic particles leading to uneven ra- 
dio emission, or the lobe volume could simply be overestimated, for 
example if their outer surfaces are more structured than a smooth 
ellipsoid and the lobes contain pockets or filaments of IGM plasma. 
In principle, it is possible to place a limit on the effects of the last 
of these possibilities from the X-ray data. If a fraction of the ap- 
parent lobe volume contains IGM gas, the expected X-ray surface 
brightness decrement will be reduced. Unfortunately the IGM sur- 
face brightness is low, only small decrements are expected, and the 
large uncertainties mean that we cannot determine a useful limit in 
this way. However, the presence of the cavity rims, probably con- 
sisting of material swept up by the expanding lobes, argues that 
the lobes have been effective in removing IGM gas from the vol- 
ume they occupy. A detection of inverse Compton emission from 
the lobes (arising from CMB photons scattered up to X-ray ener- 
gies by relativistic electrons) might allow us to place limits on the 
clumping of the relativistic plasma. Inverse Compton emission has 
been detected from the lobes of many FR II radio galaxies (Cros- 
ton et al. 2005) and from the diffuse mini-halo in the core of the 
Perseus cluster (Sanders & Fabian 2007), but is not detected in the 
lobes of bright, cluster-central FR I radio sources such as Hydra A, 
Hercules A and M 87 (Hardcastle & Croston 2010; Simionescu 
et al. 2008). Unfortunately, the inverse Compton emission from the 
lobes of NGC 4261 is expected to have only a small flux if the min- 
imum energy condition is correct, ~2 njy at 1 keV for each lobe. 
This compares with 2.6 nJy measured for the resolved X-ray jet 
(Worrall et al. 2010), but the expected 2 nJy would be spread over 



the whole area of the lobe and would be undetectable in our data. 
We are therefore unable to place limits on the filling factor, and 
must continue to assume 4>=\. 

The additional pressure required could be provided by ad- 
ditional non-radiating particles in the relativistic plasma. Since 
we require k~4, these would probably have to be introduced into 
the relativistic plasma via the entrainment and heating of thermal 
plasma by the radio jets. It has been suggested that some knots in 
the inner jets of Centaurus A may be caused by interactions be- 
tween jet plasma and the winds of high mass stars (Hardcastle et al. 
2003), and stellar mass loss from stars within the jets is an obvious 
source of entrained material. The well defined small-scale X-ray 
jets have a width of 4.7" (720 pc), and length 31.7" (4900 pc) for 
the longer western jet (Worrall et al. 2010), and we assume they are 
cylindrical. Using the mean stellar luminosity density within this 
radius (~0.11 Lb© pc -3 , based on optical modelling described in 
Section 4.2) and adopting the stellar mass loss rate from AGB stars 
M=0.0788(L B /10 10 L SQ ) M Q yr" 1 (Athey et al. 2002), we esti- 
mate the rate of mass loss into each jet to be 1 .67 x 10~ 3 Mq yr~ 1 . 
Assuming the timescale of the current AGN outburst to be our es- 
timated upper limit (75 Myr), we would thus expect a minimum of 
~ 125000 Mq of stellar material to be entrained by each jet. 

If this entrained material takes the form of a thermal plasma 
mixed with the relativistic plasma of the radio lobes, the tempera- 
ture required for it to balance the pressure of the surrounding IGM 
is >10 MeV. We note that Sanders & Fabian (2007) place an up- 
per temperature limit of 100 keV on any thermal plasma in the ra- 
dio lobes of Perseus A. For electrons, a temperature of 10 MeV is 
equivalent to a Lorentz factor 7 ~ 20. The very low density of 
such a plasma component (n e ~ 3 x 10~ 7 cm~ 3 ) means that it is 
unlikely to produce any detectable emission either directly or via 
inverse-Compton scattering. Electrons of 7 ~ 20 in the minimum 
energy magnetic field would be expected to radiate at ~ 1 kHz, be- 
low observable frequency bands. 

Limits can also be placed on the density of any thermal plasma 
component in the lobes, based on the degree of depolarization ob- 
served. The relationship between the observed degree of polariza- 
tion Tots, the the intrinsic degree of polarization Tint can be ap- 
proximated as 



Vobs. = Vint 



sm(RM\ 2 



(3) 



RMX 2 ' 

where RM is the rotation measure and A the wavelength (Govoni & 
Feretti 2004). The rotation measure of a thermal plasma intermin- 
gled with the relativistic plasma in the radio lobe is approximately 



RM = 4.05 x 10' 



B\\ dz radm 



(4) 



where n e is the electron number density of the thermal material 
in cm -3 , -By is the line-of-sight component of the magnetic field 
in Gauss, and the thickness of the region, L, is in pc (Worrall & 
Birkinshaw 2006). A depolarization of ~8 per cent, at 21.2 cm is 
reported for the lobes of NGC 4261 (Bologna et al. 1969) and while 
the contribution from beam depolarization is unknown, we can 
limit the density of any thermal component to be <2x 10~ 4 cm~ 3 . 
This is considerably greater than the density expected from en- 
trained gas from stellar winds (3.6x 10~ 7 cm~ 3 for the east lobe) 
but comparable to the IGM density. The limit thus argues only 
against the entrainment or envelopment of very large quantities of 
IGM gas into the radio lobes, but cannot rule out entrainment of 
gas from stellar winds. 
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The estimate of the mass of gas from stellar winds entrained 
by the jets ignores stars in the outskirts of the galaxy where the 
jets are broader, and the possibility of entrainment from the gas of 
the cool core or IGM. Increasing the mass entrained would lower 
the temperature of any thermal component within the lobes, but the 
energy required to heat such a component, ~5x 10 57 erg summed 
over both lobes, would not be significantly altered. This is only ~20 
per cent, of the estimated total mechanical power output of the jets, 
probably within the uncertainties on the available energy. Given the 
evidence of rims around the cavities, indicating that the lobes prob- 
ably have a high filling factor and may be mildly overpressured, it 
seems likely that the true mechanical jet power is somewhat higher 
than our estimate, and this could easily provide the energy required 
to heat any entrained material. 



5 DISCUSSION 

5.1 The cool core and AGN 

We can consider three sources of material which might fuel the 
AGN in NGC 4261: 1) a cold gas reservoir unrelated to the hot 
IGM, 2) hot gas produced within a cool corona, with no significant 
inflow of gas from the IGM to the corona, and 3) hot IGM gas 
via a cooling flow, with no separation between the cool core and 
the rest of the IGM. For AGN feedback to control cooling in the 
group, the AGN must be fuelled by material which it is capable 
of effectively heating, so that its power output directly affects the 
rate of fuel supply. This appears unlikely in case 1 and most likely 
in case 3, since we observe the interaction between the AGN and 
IGM. However, if the cool core is a galactic corona, magnetically 
separated from the surrounding IGM, then significant inflow from 
the IGM cannot occur and case 3 is ruled out. A feedback loop is 
possible in case 2 only if the AGN can effectively heat the gas of 
the corona or prevent its transport in to the central engine. 

The evidence of past tidal interaction in NGC 4261 raises the 
question of whether the current nuclear activity could be related 
to a merger. This appears unlikely. Although a nuclear disk of gas 
and dust containing (5.4±1.8)x 10 4 M of molecular and atomic 
hydrogen is observed (Ferrarese et al. 1996), cool gas has not 
been detected elsewhere in NGC 4261 (Serra & Oosterloo 2010; 
Combes et al. 2007). Stellar population modelling finds little ev- 
idence for a young stellar population in the galaxy core; the lu- 
minosity weighted age is estimated to be 12.6^'^ Gyr (Serra & 
Oosterloo 2010). It therefore seems likely that any merger was gas 
poor. 

While a large injection of cool gas into NGC 4261 is unlikely, 
it seems probable that the cool gas disk is the immediate source of 
material fuelling the AGN. Direct accretion from the hot phase can- 
not be ruled out, but does not explain the presence of the cool gas. 
Assuming a conversion efficiency of 1 per cent., accretion of the 
material in the observed gas disk would release ~10 57 erg, a factor 
of 20 less than the enthalpy of the cavities. Unless we are coinci- 
dentally observing NGC 4261 just as the AGN is about to consume 
the last of its fuel, the disk gas must be replenished by material 
cooling from the hot phase. Judging the likelihood of such a co- 
incidence is difficult. The selection of NGC 4261 as a luminous 
radio galaxy with extended lobes could bias us toward observing 
the system during the period when the lobes have had time to ex- 
pand and attain a high surface brightness, but before the jets shut 
down. Conversely, the lack of evidence of multiple AGN outbursts 
or star formation would require any gas introduced by a merger to 



be efficiently transported into the galaxy nucleus as a single unit. 
However, even if the current activity were fuelled by cold gas, the 
presence of the cool core strongly suggests that the AGN is inef- 
fective in heating its immediate neighbourhood. The presence of 
the cool core argues that interaction between the jets and gas in the 
core regions is minimal, and shock heating is weak. The total en- 
ergy released through shock heating need not be negligible, but it 
seems more likely that in this system any recent shocks have been 
driven by lobe/jet expansion rather than originating at the nucleus, 
and have occurred at radii > 10 kpc. 

The energetic conditions of the cool core (see Section 4.2) are 
similar to those described by Sun et al. (2007) for galactic coronae. 
Sun (2009) showed that for systems hosting relatively powerful ra- 
dio galaxies, coronae fall into a region of ix/iradio space sepa- 
rate from systems with larger cool cores (see their Figure 1). From 
the deprojection analysis, we estimate the 0.5-2 keV luminosity of 
the gas in the core (radius <9.8 kpc) to be 5.9x 10 40 erg s _1 . The 
1.4 GHz radio luminosity of 3C 270 is 2.4x 10 24 W Hz" 1 (Condon 
et al. 2002). This places NGC 4261 in the corona class. However, 
the cool core of NGC 4261 is relatively large compared to coro- 
nae in galaxy clusters and has a rather mild temperature difference 
with the surrounding IGM. This is expected, as corona size is de- 
pendent on the pressure of the surrounding IGM, which is lower in 
this group than in higher mass galaxy clusters. The limits on the 
suppression of conduction and efficiency of AGN heating are also 
comparatively weak. We conclude that the cool core of NGC 4261 
is probably an example of a nearby, relatively large galactic corona, 
but that we cannot be sure of the classification in such a system. 

The radio source has clearly affected the structure of the 
corona. Worrall et al. (2010) point out that the X-ray jets are sur- 
rounded by regions of low surface brightness, suggesting that ther- 
mal gas has been removed and the lobes now extend back to fill 
these regions, surrounding the jets. The X-structure observed in 
X-ray surface brightness and temperature maps, and the cool tem- 
perature of the structure, suggests that the gas of the cool core has 
been compressed and driven away from the jets by subsonic expan- 
sion of the lobes. Direct interaction with the jets appears unlikely, 
since there is no evidence of shock heating. The expansion of the 
lobes is also probably responsible for the general north-south exten- 
sion of the X-ray emission within the stellar body of NGC 4261; 
the lobes have pushed the hot IGM gas originally occupying the 
outer part of the galaxy away from the jet axis to form an annu- 
lar structure which we are viewing edge-on. It is worth noting that 
such an annulus would have a similar alignment to the disk of cool 
gas and dust which surrounds the AGN, though its scale would be 
~200 times larger. 

Since the corona shows no signs of having been heated by the 
current AGN outburst, it is unclear how nuclear activity can stop 
cooling in this region. If the AGN cannot effectively heat the core, 
it cannot form a feedback system with the IGM, and cooling and 
heating of the group gas need not be in balance. If the feedback 
relationship is to be sustained, we require a mechanism for ending 
nuclear activity. Several possibilities can be suggested: 

(i) If the AGN is fuelled by a reservoir of cool gas left over from 
a gas-rich merger, activity will cease within the next few Myr, as- 
suming a constant rate of consumption. This requires us to be ob- 
serving NGC 4261 at an unusual period of its history, and suggests 
that FR I radio galaxies are fuelled from reservoirs of cool gas, and 
is thus unsatisfying. 

(ii) A change in accretion or jet properties could alter the heat- 
ing mechanism. For example, we could be observing NGC 4261 
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during a long period in which accretion drives jets without produc- 
ing shocks in the nuclear region. At some future point, alteration of 
conditions could cause the amount of energy released in shocks to 
rise, heating the core and ending the outburst. This cannot be ruled 
out, but we note that there is evidence that the ^75 Myr period of 
activity which has formed the radio lobes has not been continuous, 
and that the jet properties have changed with time (Worrall et al. 
2010). 

(iii) The expansion of the lobes back along the line of the jets 
into the corona raises the possibility that they could reach the cen- 
tre of NGC 4261 and envelop the nucleus. This could potentially 
reduce the amount of thermal plasma able to reach the central en- 
gine. Enclosure need not lead to an immediate cessation of nuclear 
activity, since the cool gas disk would still need to be consumed, 
and the enclosed region could also contain a high-pressure core of 
hot gas. However, it is unclear whether the pressure in the lobes will 
be sufficient to push aside the thermal gas at very small radii. As 
the lobes expand outward, moving down the IGM pressure gradi- 
ent, we might expect their internal pressure to decline. In this case, 
they are unlikely to be able to move aside the high-pressure ma- 
terial in the inner part of the corona. It is also unknown whether 
flows of cooling material might be able to penetrate or push aside 
the enclosure formed by the radio lobes, transporting gas in to the 
nucleus. 

(iv) Disturbance of the structure of the corona by the radio 
source could affect the AGN in a number of ways. The lobes have 
already compressed the gas in the core, increasing its density and 
cooling rate. This could lead to an increase in the supply of gas 
to the AGN and potentially an increased heating rate. Disturbance 
of the corona structure may also affect the conduction rate. The 
surface of the corona has already been increased above that of a 
sphere. If enough of this surface is in contact with the IGM (rather 
than the radio lobes) the increased surface area will lead to an in- 
crease in the conduction rate. Gas motions introduced as the lobes 
push the core gas aside may also alter the magnetic field structure, 
potentially increasing conduction if field lines are straightened by 
expansion, or reducing conduction if the motions cause further tan- 
gling of the field. If the suppression factor is significantly reduced, 
the corona could potentially be heated fairly rapidly, reducing the 
fuel supply to the AGN. 

While we cannot determine whether any of these processes 
(or some other) will be effective in terminating the AGN outburst 
in NGC 4261, it is likely that physical changes in the structure of 
the core over a long timescale are an important factor. If this is 
true of other corona systems, it would suggest that whereas AGN 
in large cool cores may rapidly reheat the cooling gas which fuels 
them, AGN in coronae can remain active for longer periods without 
disrupting their fuel supply. Coronae therefore seem likely to fuel 
long-term nuclear activity, rather than short outbursts, and AGN 
in coronae may spend a larger fraction of their time active than 
AGN in large cool cores. This could explain why such systems lack 
large-scale cool cores; they inject more energy into the IGM, and 
heat it to a point where the cooling which occurs while the AGN 
is quiescent is insufficient for a cool core to develop before the jets 
restart. Studies of other nearby corona systems are clearly required 
to determine whether such important differences in duty cycle and 
energy output are real. Dunn et al. (2010) examined 18 nearby gi- 
ant elliptical galaxies and found that 17 host nuclear radio sources, 
potentially lending weight to the idea that AGN in galactic coronae 
are able to remain active for long periods. However, a number of 
galaxies in the sample occupy the centres of groups or clusters, and 



have cool cores considerably larger than that of NGC 4261, so only 
a fraction of this sample are likely to be corona systems. 



5.2 Structure of the radio source 

The structure of the radio lobes, cavities and rims raises several 
questions. Considerably more X-ray emission is seen in the east- 
ern lobe and its rim appears, at least in the XMM-Newton images, 
to be thicker. It is also complete, unlike the western rim, which is 
undetected in one quadrant. Two possible explanations involve the 
data quality; the XMM-Newton data may simply be insufficiently 
deep to clearly detect the cavity rims and determine their shape, and 
our point source subtraction may also be ineffective. Prior studies 
of the X-ray point source population have shown a particular con- 
centration of sources to the north and northeast of the AGN (Gior- 
dano et al. 2005), but also to the south and west. The poorer ability 
of XMM-Newton to resolve these sources compared to Chandra 
might produce increased noise in the northeast quadrant, leading 
us to overestimate the thickness of the rim in this area. 

Other factors which could affect rim strength are the struc- 
ture of the IGM into which the radio lobes have expanded, and the 
structure of the lobes themselves. Surface brightness modelling of 
the group with the ROSAT PSPC found the IGM to be centred to the 
east of NGC 4261, with the emission more extended to the south 
and east. Davis et al. (1995) found the centroid of diffuse emission 
to be R.A. 12 h 19 m 27 s , Dec. +05°49'58" (J2000), ~1' (9.2 kpc) 
ENE of the nucleus (see also their fig. 1). If this extension is real, 
it suggests that NGC 4261 is slightly offset from the centre of the 
group X-ray halo. The eastern lobe would, in this scenario, have 
expanded into a denser environment than its western counterpart, 
and has therefore accumulated more gas in its rim. 

The morphologies of the two radio lobes are also slightly dif- 
ferent. The western lobe is broader but less extended, and has dis- 
tinct bulges on its sides. If the extension of the lobe tip is driven 
by the jet, the shorter length of the western lobe could indicate that 
the western jet has encountered greater resistance to its expansion. 
However, we see no evidence of any feature which could be re- 
sponsible for this resistance, such as a region of denser intra-group 
gas. The greater sideways expansion of the western lobe suggests 
that external pressure is lower on this side of the core, as expected 
if the galaxy is offset from the group centre. There is also some 
indication that the eastern radio jet bends inside the radio lobe (see 
the contour map in Figure la and fig. 1 of Worrall et al. 2010), 
while the western jet remains on approximately the same axis as 
the small-scale jets visible in the X-ray. This could perhaps affect 
the apparent structure of the eastern lobe, since we are viewing a 
projected image of a three dimensional structure. If the jet axis is no 
longer in the plane of the sky, the shell of gas around the expanding 
lobe tip might appear as a broader X-ray bright region rather than 
a relatively thin shell. This would probably require a fairly strong 
bend in the jet which, while quite possible, has yet to be confirmed 
from the available observational data. 



6 CONCLUSIONS 

We have analysed deep Chandra and XMM-Newton observations 
of the nearby group-central elliptical NGC 4261. The X-ray ob- 
servations reveal a large degree of structure in the gas in which the 
galaxy is embedded. The AGN jets have inflated two large lobes, 
excavating cavities in the IGM and building up rims of compressed 
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hot gas which almost entirely enclose the lobes. The western cav- 
ity has an apparent break in its rim, while the eastern lobe has a 
higher X-ray luminosity over most of its area. The enthalpy of the 
lobes is large, ~2.4x 10 58 erg, and the likely timescale of the out- 
burst, 5J75 Myr, suggests that the mechanical power of the jets, 
^10 43 erg s _1 , greatly exceeds their current radio power. The me- 
chanical power of the AGN also exceeds the X-ray luminosity of 
the gas in the region the lobes inhabit, the central 45 kpc of the 
group. The available X-ray data are insufficient to determine con- 
clusively whether the cavity rims contain shocks, but we can limit 
the maximum expansion velocity of the lobes to be M < 1.05, 
and their minimum expansion velocity to be a large fraction of the 
sound speed. These results suggest that the radio source may still 
be mildly overpressured relative to it environment, but that it is only 
inefficiently heating the IGM. 

Spectral mapping and radial spectral analysis confirms that 
NGC 4261 has a cool core with radius ~9.8 kpc and a typical 
temperature of ~0.6 keV, compared to a peak IGM temperature of 
~1.6 keV. The core is not spherically symmetric, but contains an 
X-shaped structure, with the arms of the X extending to either side 
of the radio jets, the northeast and southwest arms being most ex- 
tended. This structure was commented on by Worrall et al. (2010), 
who noted that the radio jets appear to have expelled gas in wedge- 
shaped regions to their north and south, probably corresponding 
to conical volumes. There is also a larger, bar-shaped structure 
extending north-south across the centre of NGC 4261, which is 
poorly correlated with the galaxy major axis but perpendicular to 
the jet axis. This is probably an annulus of gas viewed edge on, 
pushed into its current configuration by the expansion of the radio 
lobes. It is clear from these structures that the jets have played an 
important role in determining the structure and properties of the gas 
in and around NGC 4261. 

The properties of the cool core suggest that it should be classi- 
fied as a galactic corona. The relative radio power of the galaxy and 
luminosity of the cool core fall in the range occupied by corona sys- 
tems. Heating of the corona by the AGN and conduction must have 
relatively low efficiencies. We limit the AGN heating efficiency 
to be <3.5 per cent, and find that conduction into the cool core 
from the IGM is suppressed by a factor >5. However these values 
probably overestimate the effectiveness of the heating mechanisms 
since they assume a timescale based on the length of the AGN out- 
burst rather than the lifetime of the temperature gradient. Super- 
novae within the corona release ~6x 10 40 erg s _1 (compared to the 
0.3-7.0 keV X-ray gas luminosity of ~10 41 erg s _1 ), but since the 
stellar distribution is less centrally peaked than the gas density, they 
are unlikely to be able to prevent cooling in the inner core. Material 
lost from stars can balance the cooling rate of the corona, suggest- 
ing that cooling could fuel the AGN over timescales of several Gyr. 

Since the AGN appears unable to effectively heat the corona, 
we have considered alternative mechanisms by which cooling could 
be disrupted or stopped. We conclude that changes in the corona 
structure, caused by the intrusion of the radio lobes into the cool- 
ing region, may be able to alter the cooling rate by compressing 
the gas, isolating the nucleus, or increasing the conduction rate. In 
any case, these results suggest that corona systems can potentially 
sustain longer periods of nuclear activity compared to their coun- 
terparts in large cool cores. Since the powers of radio galaxies with 
coronae and large cool cores are similar, it seems likely that sources 
in coronae will, on average, release more energy into the IGM over 
long timescales. Given the large numbers of coronae already de- 
tected in group- and cluster-central galaxies (e.g., Sun 2009), this 



may be an important consideration in determining the gas temper- 
ature structure of the population of galaxy systems. 

Comparing the pressures of the IGM gas and the relativis- 
tic plasma of the lobes, we find the apparent lobe pressure to be 
a factor of 5 lower than that of the surrounding halo. This indi- 
cates that either additional non-radiating particles are present in 
the lobe plasma, that the filling factor of the lobes is low, or that 
the magnetic field departs from the minimum-energy solution. The 
currently available data are insufficient to place useful limits on the 
filling factor of the lobes. Thermal gas entrained from stars within 
the radio jets, or from the IGM, could potentially provide the addi- 
tional pressure without radiating in a detectable waveband or pro- 
ducing excessive polarization effects. The energy required to heat 
this gas to its presumed temperature is ~20 per cent, of the total 
enthalpy of the radio lobes. 
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